Topological invariants for interacting topological insulators with inversion symmetry 
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For interacting Z2 topological insulators with inversion symmetry, we propose a simple topological 
invariant expressed in terms of the parity eigenvalues of the interacting Green's function at time- 
reversal invariant momenta and zero frequency. We derive this result from our previous formula 
involving the integral over the frequency- moment a space. This formula greatly simplifies the explicit 
calculation of Z2 topological invariants in inversion symmetric insulators with strong interactions. 
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I. INTRODUCTION 

Topological insulators are quantum states of matter 
with insulating bulk and stable metallic surface^"— . The 
origin of their interesting physics is the bulk topology in 
momentum space, which can generally be characterized 
by topological invariants, the simplest examples among 
which are the TKNN invariant^ and the Z2 topologi- 
cal invariants''"^^. The topological band invariants are 
defined in terms of the Bloch states, therefore, strictly 
speaking they can apply only to non-interacting systems, 
although one may extrapolate them to weakly interact- 
ing systems. On the other hand, the axion angle in the 
electromagnetic response can provide a general defini- 
tion of an interacting topological insulator—. In recent 
years, after the prediction and discovery of weakly inter- 
acting topological insulators, topological insulators with 
strong interaction have been among the main topics in 
the field—"—. It is therefore urgent to search for topo- 
logical invariants whose applicability is beyond the the 
non-interacting limit. One of the promising directions 
is the Green's function approach^!, which evaluates the 
axion angle and remains valid in the presence of interac- 
tion, provided that we use the full single particle Green's 
function. Much recent interest has been focused on this 
approach^^"--. 

Motivated by the topological field theory approach and 
the idea of dimensional reduction^i^, we previously pro- 
posed a Z2 topological invariant in terms of Green's func- 
tion in the extended frequency-momentum space^i: 
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in which the Greek letters are four-dimensional 
momentum- frequency indices including po = i'^i and u 
is an extension parameter between the physical Green 
functions at u = and a trivial constant reference func- 
tion aX u — \. For different extrapolations the value of 
P3 obtained may differ by integer, but the fractional part 
of P3 mod 1 is independent from the extrapolation and 
is a physical topological invariant. This topoological in- 
variant has the form of a Wess-Zumino-Witten (WZW) 



term, which is not easy to implement in practical calcu- 
lation. In the present paper, we propose much simpler 
form of topological invariants [Eq.JS]) and Eq. pH]) ] for 
interacting topological insulators with inversion symme- 
try. In this new topological invariant. Green's function 
at eight time reversal invariant momenta (TRIM) Ti is 
evaluated and there is no integral involved, which greatly 
simplifies the calculation. This topological invariant can 
be regarded as a highly nontrivial generalization of the 
Fu-Kane invariant^^ for non-interacting system. From 
the mathematical perspective, the logic leading to our 
formula is analogous to that of refi^. 



II. TOPOLOGICAL INVARIANT FOR 
THREE-DIMENSIONAL TOPOLOGICAL 
INSULATORS VIA GREEN'S FUNCTION 
EIGENVECTORS AT ZERO FREQUENCY 

In this paper we focus on three-dimensional time- 
reversal invariant insulators with inversion symmetry. 
The time reversal operator is an anti-unitary operator 
defined by T = TK, where K is the complex conjuga- 
tion operator and T is a matrix satisfying T*T ~ — 1. 
The inversion transformation matrix P satisfies ~ 1, 
therefore, its eigenvalues are ±1. The Green's function 
is a A'^ X matrix, where N is an even integer due to 
time reversal symmetry. It is convenient for our purpose 
to use the imaginary frequency Green's function, e.g. a 
free fermion system with Hamiltonian H ~ X^fc Cfc^(fc)cfc 
has Green's function G{iu>,k) = l/[iuj — h{k)]. Gener- 
ally, we have the Schwinger-Dyson equation G^^{iuj, k) — 
GQ^{iuj,k) — S(ia;,fc) between the full single particle 
Green's function G, the non-interacting Green's func- 
tion Go, and the self energy E. A general G{iuj,k) G 
GL(iV, G) can be diagonalized to give the eigenvalues and 
eigenvectors. Explicitly, we have 

G~^ [iuj , k)\a{iu) , k)) — iia{i^,k)\a{iuj,k)) (2) 

where fiaiii-^, k) = aa{iuj, k) -\- iba{iio, k), with a^, ba be- 
ing real number. The matrix G{iuj, k) has the same 
eigenvectors but eigenvalues iJL~^{iu),k). Because of in- 
version symmetry. Green's function satisfies PG{iuj, k) — 
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G{iuj, -k)P. At eight TRIMs k = fj, where is equiv- 
alent to — Fi, the Green's function satisfies PG{iuj,Ti) = 
G{iuj,ri)P, therefore, \a) can be chosen as the simulta- 
neous eigenvectors of G and P, namely that in addition 
to eq.([2]), I a) also satisfies the equation 

P\a{iuj,fi)) ^ r]a\a{iuj,fi)) (3) 

where rja is the parity eigenvalue. Using the Lehmann 
representation, we can obtain the equation 

G^{-iuj,k) = G{iuj,k) (4) 

from which it follows that G(0, F^) is a Hermitian matrix 
and thus has real eigenvalues, i.e. ^^(O,?^) is always lo- 
cated on the real axis of the /Xq, = Ua+iba complex plane. 
If fia{0,^i) is located on the right half of the real axis, 
namely iJ,a{0,^i) > 0, we call \a{0,Ti)) a "Right-zero 
(R-zero)" . Similarly we can define "L-zero" . One of the 
central results of this paper is the topological invariant 
defined as 

R— zero 

where we have used the convention that (—1)^/^ = i 
for the square root. The appearance of square root 
can be interpreted as follows. From the equations 
r-iG'^(iw,fi)T = G{iuj,ft)^^ (the superscript is 
matrix transposition) and Eq.Q, we can show that 

G-\-iio, f,)T\a{iLJ, f ,))* = ^il{^to, f,))T|a(^c^, f,))*(6) 

Therefore, R-zeros |a(0, fi)) and r|a(0,fi))* always 
form a Kramers pair. The square root in eq.© amounts 
to requiring that each Kramers pair with rja — con- 
tributes = —1. We will provide a derivation of eq.® 
in the following sections, and also present an equivalent 
formula in eq.©. 

In the non-interacting limit, the parity of Green's func- 
tion eigenvectors is reduced to the parity of Bloch states, 
and it can thus be shown that our formula in eq.® is 
reduced to the Fu-Kane formula3^. 



III. AN EQUIVALENT FORMULA AND THE 
TOPOLOGICAL INVARIANCE 

Let us begin with some notations. As uj flows from — oo 
to -|-oo, /Xct(ia;,Fi) sweeps a curve on the complex plane. 
There are two special points and oo on this complex 
plane, which can be compactified by adding a point oo 
at infinity. We assume that G{iui, F^) G GL(7V, C), there- 
fore, /ia 7^ Ooroo. Because ^aii^, k) ^ iuj as uj ±oo, 
each fj,a curve must cross the real axis by an odd number 
of times. All /Xq, curves can be classified into two topo- 
logical classes: "right fia curves" (R-curves) and "left 
curves" (L-curves). The R/L-curves, by our definition, 




(c) 



FIG. 1: Topological invariance of R/L-curves on the fia com- 
plex plane, (a) A and B is L-zero and R-zero respectively, 
and the corresponding curve is L/R-curve respectively, (b) 
Compactified complex plane, on which a R-curve I is shown. 
Here the only crossing point of I with real axis is D. The point 
oo is not reached by although I approaches cxd as cj — >■ ±oo. 
(c) The R-curve h can be deformed into another R-curve h 
following the path indicated by the arrow, however, it can- 
not be deformed into a L-curve if and oo are avoided, i.e. 
topological class is preserved during smooth deformations. 

cross the right/left half real axis at iuj = 0[see Fig[T^]. It 
is also worth noting that the fia curve cross the real axis 
only at zo; = 0, which follows from the Lehmann repre- 
sentation. Due to Eq.Q, each R/L-curve is symmetric 
according to the real axis [see Fig.([T^)]. Furthermore, 
from Eq.® we see that each R/L-curve has a Kramers 
partner which coincides with it. It should be empha- 
sized that the distinction between R-curves and L-curves 
is topological because an R-curve cannot be smoothly 
deformed to a L-curve without crossing the two singular 
points and oo, nor vice versa [see Fig. ([!])]. 

Now we can define a new topological invariant in terms 
of R/L-curves as 

(-1)^= n ^^'-±1 (7) 

R— curve 

which is equivalent to eq.(IS]), as we will see. The parity 
formulas in eq.® and eq.© are the central results of this 
paper. The topological invariance of eq.([7]) follow from 
the topological invariance of R/L-curves and the fact that 
the parity rja, which takes only two discrete values ±1, 
cannot changes abruptly during a smooth deformation. 
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It is clear that eq.(l5]) is more efficient in practical calcu- 
lation, while eq.(I71) provides a deeper perspective, for in- 
stance, we note that the definition of "R/L-curve" can be 
modified as follows: We can arbitrarily choose a straight 
line (other than the imaginary axis ) through the singu- 
larity 0, and define "R/L-curve" according to the loca- 
tions where fia curves cross this line. It can be easily 
appreciated that the topological class of an R/L-curve is 
independent of the choices of this straight line. Eq.([7]) 
reveals such degree of freedom in the definition of R/L- 
curve. 

We have assumed that Green's function has no singu- 
larity, which is a mild assumption for an insulator. How- 
ever, if we want to study phase transitions, singularities 
must be considered. There are two types of singularity 
of Green's function, one of which is the crossing of /Zq 
curves with 0, the other is crossing with oo. Topological 
phase transitions must be accompanied by these singu- 
larities. The transitions of first type (singularity at 0) 
include all topological phase transitions in free fermion 
systems, while those of the second type (singularity at oo) 
can happen when the interaction is sufficiently strong. 

IV. DERIVATION OF THE NEW FORMULA 
FROM THE PREVIOUS P3 FORMULA 

In this section we provide a derivation of the parity 
formula in eq.(IS|) and eq.©. It was shown in rcf.^^ that 
the Z2 topological invariant of interacting topological in- 
sulators can be expressed in terms of Green's function as 
in Eq. ([T|). It can be expressed more compactly as 

2P3 = WiG)\R^T^ 

1 r"^ 

- 4807r3 ^ 

xGd^G^^GdpG-^GdaG-^GdrG-^] 
= integer (8) 

where Vl^(G)|flxT'i is the "winding number" of the map 
from frequency-momenta space RxT^ to GL{N, G). We 
have used the conventions that p — (po,pi,p2,P3,P4), 
Po = iuj is the frequency, {pi,P2iP3) are momenta, and 
P4 is the WZW parameter. (po,Pi,P2,P3, 0) expands the 
physical frequency- moment a space, and (pojPijP2iP3j tt) 
expands the frequency-momenta space of the trivial ref- 
erence system. We have assumed that the eigenval- 
ues of Green's function is everywhere smooth and non- 
vanishing. The topological classes of Z2 insulators are de- 
fined by 2P3 mod 2. For inversion symmetric system, we 
can also add the constraint that PG{po,Pi,P2TP3,Pi) = 
G{pq, -pi, -P2, -P3, -Pa)P- 

With this preparation, we would like to obtain the 
equation 2P3 = A mod 2, in which 2P3 is defined by 
eq.® and A is defined by cq.© or eq.©. The basic 
idea is outlined as follows. Both 2P3 and A are topologi- 
cally invariant, therefore, as we deform G{p) a.s a. N x N 
matrix, 2P3 and A can only change abruptly when G{p) 



encounters singularities. We can compute their changes, 
and we will find that they are equal mod 2. Now we 
only need one example in which 2P3 — A mod 2 to prove 
this identity in general cases, but we already know that 
2P3 — A mod 2 in the non-interacting limit^*^'^^, there- 
fore, 2P3 = A mod 2 is generally correct. 

Now we present our calculations following the above 
outline. When we talk about deformation of G{p), 
the precise meaning is that G = G{p,X), where A 
parametrizes the deformation. Without losing generality, 
we assume the singularity occurs at A = 0. We consider 
a singularity located at p — (0, 0, 0, 0, 0) for the following 
reasons (for the spatial components, pi = 0; z = 1,2,3 
means the TRIM). The spatial components {pi,P2,P3) 
are assumed to be one of the eight TRIMs, because 
singularities away from TRIMs must occur pairwise at 
(pi,P2,P3,P4) and (-pi, ~p2, -P3, -P4) due to inversion 
symmetry, which can thus only change the value of 
W{G)\fixT'^ by an even integer. Similarly, the frequency 
Po — because of the equation G^ {—pojTi) ~ G{po,ri), 
which requires that the singularities occur pairwise at po 
and — Po- 

First we consider the singularity defined by 
det[G-^p,\ = 0)] = 0. We can expand G-^p,X) 
in a neighborhood keeping only terms linear in p. In 
the cases when the first order terms vanish and higher 
order terms dominate the expansion, we can always 
perturb the singularity to split a single singularity into 
several linear ones. At a generic linear singularity, four 
eigenvalues of G~^{p,X = 0) vanish, as is required by 
time-reversal and inversion symmetries. Therefore it is 
sufficient to write only the singular part of the Green's 
function, which has the following form: 

4 

G"^(p, A) = Uo{po,\) + '^Ua{pt)Ta +U5(P0,A)P (9) 

where uo(Po,A) = CqoPq + G05X, ^5(^0, A) = C50P0 + 
C55A, and Ua = Y^*=iGazPi{a = 1,2,3,4). Generically 
we have det G ^ 0. The matrices appearing here are 
P = 1 (g)T3, Fq = cTa ®Ti{a = 1, 2, 3), F4 = 1 (g)T2, whcrc 
both a and t are Pauli matrices. Time reversal matrix is 
defined as T = ia2 1 in our notation. The form of this 
expansion is dictated by PG{po,Pi) — G{po, —pi)P and 
T^^G^{po,Pi)T = G{pa,-pi), e.g. inversion symmetry 
requires that Fq {a = 1, 2, 3, 4) contain factor ti or T2 but 
not T3 . We note that the expansion in eq. (|9]) is analogous 
to that in refi^i^, in which an expansion of Hamiltonian 
according to symmetries was used. 

We now proceed to calculate the change W{G, X = 
0+)|flxT4 - W{G,X = 0-)UxT4. To this end, we de- 
fine a five dimensional sphere by X]^=o b^P + = 
r^, whose radius r is sufficiently small, then we have 
a mapping / from to the manifold parameterized 
by {(uo. Ml, U2, 1*3, U4, U5)}. From a simple geometri- 
cal picture[see Figl2], we have W{G,X = 0+)|flxT-i — 
W{G,X = 0-)|^xT4 = W{G)\s5. For the standard 
mapping (po,Pi,P2,P3,P4, A) (uq, ui, U2, U3, "4, W5) = 
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ally reduced to three dimensions, the robust topological 
invariant is exactly the Z2 part. 

V. INTERACTING TOPOLOGICAL 
INVARIANTS FOR TWO-DIMENSIONAL 
TOPOLOGICAL INSULATOR WITH INVERSION 
SYMMETRY 

For two dimensional time reversal invariant and in- 
version symmetric topological insulators, similar formu- 
las like eq.® and cq.© can also be defined, with the 
only difference that there are now four instead of eight 
Ti points. Explicitly, we have 



FIG. 2: The deformation of integral contours. The p axis 
refers top = (po,Pi,P2,P2,P4). The two contours in (a) corre- 
spond to A > and A < respectively. The opposite orienta- 
tions come from the negative sign in —W{G, A = 0~). We can 
deform these two contours to those illustrated in (b), which 
does not change W(G) due to topological invariance. In the 
calculation of the jump of W{G), the "dashed lines" (which 
are in fact 5d manifolds) in (b) cancel each other, and we 
are left with the "solid circle", which is the five dimensional 
sphere S^. 



(-1)^^°= n ^'J'- n ^'J' (10) 

R— zero R— curve 

whose topological invariance and derivation follows the 
same logic as the previous section, which we shall not 
repeat. 



VI. CONCLUSIONS 



^{po,Pi,P2,P3,P4, X), where i? is a nonzero constant, 
we can obtain W{G)\s5 1 by strightfoward calcula- 
tions. Because of the linear relation between Ua {a = 
0,1,2,3,4,5) and pQ,pi,p2,P'i,PA,X, we can obtain the 
general result VF(G)|55 — ±1. 

At the same time, we can see that the R/L-zeros have 
tta = Re(Co5 + C5577q)A, rja — ±1, both of which are 
two-fold degenerate. We look at the two eigenvectors 
with ?7c( = — 1 because only those contribute nontriv- 
ially to A. When A changes sign, the two R/L-zeros 
with ?7a = — 1 will change their topological classes, i.e. 
from being R-zeros to being L-zeros or vice versa. Cor- 
respondingly, two R/L-curves with r]a = —1 change their 
topological classes. This will generate a factor = ~1 
m (-1)'^. Therefore, we see that the changes of 2P3 and 
A match each other. In the above analysis we consider 
an singularity with vanishing det [G~^(p)], but we can 
also consider singularities with vanishing det [G(p)]. The 
linear expansion around the singularity is the same as 
that in vanishing det [G~^(p)] cases. In summary, we 
show that for every linear singularity, the changes of 2P3 
and A are always equal mod 2. Since all other types of 
singularity can generically be perturbed to linear ones, 
we have proved our statement that 2P3 = A mod 2. 

By an analogous calculation, one can show that a 
four (spatial) dimensional time reversal invariant insu- 
lators with inversion symmetry is topologically nontriv- 

ial if riR-zero^"^^ ~ ^^^^ parity approach, only 

the Z2 part of Z classification can be detected. When 
the four-dimensional topological insulator is dimension- 



In this paper we proposed a simple form of the Z2 
topological invariant for interacting topological insula- 
tors with inversion symmetry. These formulae are well 
defined provided that the Green's function has no singu- 
larities. In these cases, our formulae provide a Z2 clas- 
sification of interacting topological insulators. It should 
be mentioned that this is not a complete classification 
of time reversal invariant insulators because those states 
with ground state degeneracy are not included, as men- 
tioned in Refj^. However, most of the realistic models 
and materials with inversion symmetry are within the 
description presented in this paper, thus the new formu- 
lae can greatly simplify the evaluation of Z2 topological 
invariants. 

In the inversion-asymmetric case, if one can tune phys- 
ical parameters so that the system smoothly evolves into 
a inversion-symmetric one, i.e. the Green's function is 
nonsingular during the evolution, then the Z2 class of the 
original system can be inferred from that of the inversion- 
symmetric system, which can be calculated using our 
parity formula. This parity formula can be calculated 
explicitly using various numerical methods. It is also 
worth mentioning that our formula is also applicable to 
the so-called "inversion symmetric topological insulator" , 
i.e. topological insulators protected only by inversion 
symmetry ^-^"^^ . 
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